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Abstract 

In the multiple monomers per site (MMS) model, polymeric chains are represented by walks on 
a lattice which may visit each site up to K times. We have solved the unrestricted version of this 
model, where immediate reversals of the walks are allowed (RA) for K = 3 on a Bethe lattice with 
arbitrary coordination number in the grand-canonical formalism. We found transitions between a 
non-polymerized and two polymerized phases, which may be continuous or discontinuous. In the 
canonical situation, the transitions between the extended and the collapsed polymeric phases are 
always continuous. The transition line is partly composed by tricritical points and partially by 
critical endpoints, both lines meeting at a multicritical point. In the subspace of the parameter 
space where the model is related to SASAWs (self-attracting self-avoiding walks), the collapse 
transition is tricritical. We discuss the relation of our results with simulations and previous Bethe 
and Husimi lattice calculations for the MMS model found in the literature. 
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I. INTRODUCTION 



The thermodynamic behavior of polymers, both in solution or in a melt, may be studied 



using continuum or lattice models 



1]. In particular, linear polymers in lattice models are 



usually described by self- and mutually avoiding walks on the lattice (SAW's). The excluded 



a. 



volume interactions, are essential to reproduce the correct scaling behavior of the system 
This constraint also adds considerable difficulties to the study of the models, when compared 
to the case of random walks, where much statistical results are known analytically [3]. As 
an example of the effect of the self-avoidance constraint on the asymptotic properties of 
a single walk on a lattice, we may recall that, while the size of the region occupied by a 
random walk with I steps on a lattice, measured by the end-to-end distance or the radius of 
gyration, grows as £ 1//2 in the limit I — >• oo, for lattices of dimension below 4, the asymptotic 
behavior for SAW's is described by an exponent which is larger than 1/2, and thus the size 
of the region occupied by the walk on these lattices grows faster with the number of steps 
of the walk when excluded volume interactions are present. For two-dimensional lattices 
this exponent is known to be equal to 3/4 [4J. We notice that since this exponent is larger 
than 1/d = 1/2 in this case, the density of monomers vanishes in the region occupied by the 
polymer. 

The basic property which describes the behavior of a single self-avoiding walks on a 
lattice is the number of walks with £ steps, starting from the origin of the lattice. We may 
consider the walks to be chains, so that the steps are bonds which link successive monomers 
of the polymeric chain. The number of monomers of a chain, which may be called its 
molecular weight M, is the number of lattice sites visited by the SAW, so that M — i + 1. 
If we wish to study a single chain in the grand-canonical ensemble, where the number of 
monomers fluctuates, we associate a fugacity z = exp (/?//) to each monomer in the chain, 
where = (fc^T) -1 and /i is the chemical potential of a monomer. The grand-canonical 
partition function will then be given by: 

Y(z) = J2Cmz M , (1) 

M 

where Cm is the number of configurations of a chain with M monomers (M — 1 steps). Alter- 
natively, this partition function may be viewed as the generating function for the numbers of 
chain configurations Cm- So far, all allowed configurations are associated to the same energy, 
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and thus the model is athermal. The model defined in this way displays a phase transition, 
a non-polymerized phase is stable at low values of the fugacity z, and for fugacities above a 
critical value z c a polymerized phase is stable, with a positive density of monomers placed 
on the lattice. At the critical fugacity, the density of the polymerized phase vanishes, so 
that the polymerization transition is continuous. The critical value of the fugacity is related 
to the asymptotic behavior of the numbers of SAW's Cm in the large M limit. There is 
good numerical evidence that Cm ~ M 1 ~ l q^ 4 , where the effective coordination number q e 
is smaller than the coordination number of the lattice, and the critical exponent 7 is equal 

n 

to 4/3 in two dimensions, 7/6 in three dimensions and 1 in four dimensions or above [5[. 
The effective coordination number is the inverse of the critical fugacity q e = l/z c , and it is 
easy to show that the grand-canonical partition function Eq. (TjQ) is singular at this value of 
the fugacity, its asymptotic behavior being given by Y(z) ~ A(l — q e z)2. In the canonical 
ensemble, the system is critical in the thermodynamic limit M — > 00 [6j. The recognition 
that the contributions to the high-temperature series expansion of the n-vector model of 
magnetism in the limit n —> reduce to SAW's on the lattice (7'j has allowed the application 
of renormalization group methods to the polymer transition, linking this problem to the 
much studied ferromagnetic transition in the n-vector model. 

The athermal polymerization model may be generalized by including attractive interac- 
tions between monomers located on first neighbor sites and which are not connected by a 
polymer bond. This model of self-avoiding self-attracting walks (SASAW's) is usually used 
as an effective model to study the behavior of a polymer in a poor solution, the attrac- 
tive interaction mimics the energetically unfavorable contact between solvent molecules and 
polymeric monomers Now, besides the monomer fugacity, an additional parameter is 
present in the model, the Boltzmann weight u = exp(/3e), where — e is the energy associ- 
ated to each monomer-monomer interaction, and there will be a competition between the 
repulsive excluded volume interactions and the attractive interactions. The model reduces 
to the previous one for 00 — 1, and as u is increased the continuous polymerization transition 
happens at lower values of the fugacity z, becoming discontinuous if uj exceeds a value ujtc- 
Thus, a tricritical point is found in the phase diagram of the model at (ztcUtc), as ma y 
be seen in the schematic diagram shown in Fig. [TJ In the canonical situation, the system 
is always on the border between the non-polymerized and the polymerized phase [6|. At 
high temperatures (low values of uj), the polymerized phase is indistinguishable from the 
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non-polymerized phase, and thus has vanishing monomer density. This phase is sometimes 



called coi 



phase in the polymer physics literature. Below the tricritical temperature (called 



6 point the chain is collapsed and the polymerized phase has nonzero density (globule 
phase). Again it is possible to map the SASAW's model on a generalized ferromagnetic 
n- vector model [8|]. In two dimensions, the tricritical point of the model has been studied in 
detail using transfer matrix techniques jsj, and the tricritical value of the exponent which 
characterizes the scaling behavior of the radius of gyration is = 4/7 10|. 
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FIG. 1: (Color online) Schematic phase diagram of the model of self-attracting self-avoiding walks 
(SASAW's) on a lattice. The continuous polymerization transition is represented by a full line 
(black on line) and the discontinuous transition is located at the dashed line (red on line). Both 
transition lines are separated by a tricritical point, represented by the circle (green on line). 



More recently, an alternative model has been proposed by Krawczyk et al [ll| for the 
collapse transition of polymers. In this model, which we may call MMS (multiple monomers 
per site) model and is a generalization of the Domb- Joyce model 12|, up to K monomers 
may occupy the same site of the lattice. The canonical version of the model, with chains of 
fixed (large) number of monomers M, was studied for K = 3 using computer simulations on 
the square and cubic lattices. Besides the case with no additional restrictions, which was 
named RA (immediate reversals allowed) model by the authors, a more restrictive model, 
where the chain is not allowed to return to the lattice site it occupied two steps ago, (RF 
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model), was also studied. Collapse transitions were found only for the RF model on the 
cubic lattice, indicating that, at least for this lattice, the restrictions seem to be essential for 
the existence of these transition. The weight of a site with two and three monomers in the 
model is oji and U2, respectively, and in the two-dimensional parameter space of the model 
defined by the variables /?i = Inui and 02 — hi 1^2 > it seems that the collapsed polymerized 
phase (globule) is separated from the regular polymerized extended phase (coil) by lines of 
continuous and discontinuous transitions, both transition lines meeting at a tricritical point. 
We note that, in the SASAW's model discussed above, the extended-collapsed transition in 
the canonical situation is continuous and of tricritical nature. One point which needs to be 
understood is the apparent absence of transitions in both models on the square lattice and 
in the RA model on the cubic lattice. 

Motivated by the questions above, the grand-canonical version of the MMS model was 
studied on Bethe and Husimi lattices. Initially, both versions of the model (RA and RF) with 
K = 2 were solved on a Bethe lattice with general coordination number q fl3j ]. Although 
these initial calculation for the models resulted in phase diagrams with some qualitative 
differences when compared to the usual behavior of SASAW's, a revision using a different 
and better fundamented procedure to find the coexistence loci resulted in diagrams which 



are similar to the ones for SASAW's, in both cases (RA and RF 



with continuous transitions 



between the polymerized phases in the canonical formalism [14j . The solution of the K = 2 
RF model on the Husimi lattice 15|] lead to a phase diagram similar to the one found for 
the same model on the Bethe lattice. A natural interpretation of this model is to consider 
that monomers on the same lattice attract each other, so that the statistical weights of sites 
with one and two monomers will be u± = z and Co>2 = ujz 2 , where z is the fugacity of a 
monomer and u = exp(— /3e) is the Boltzmann factor associated to the attractive interaction 
energy e between monomers on the same site. While the tricritical point in the Bethe lattice 
solution of the model corresponds to utc = 1, the solution on the Husimi lattice shows 
the tricritical point located at to w 1.09, in the region of attractive monomer- monomer 
interactions, as expected. More recently, the RF model for K = 3 was solved on the Bethe 
lattice in the grand-canonical ensemble {(|. In the two-dimensional subspace of the three- 
dimensional parameter space used in the grand-canonical calculations, which corresponds 
to the canonical phase diagram, again the transitions between the polymerized phases are 
always continuous. Two transition lines, one composed by tricritical points and the other 
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by critical endpoints, meet at a multicritical point, not far from the region in the parameter 
space where the tricritical point was found in the original simulations of the K = 3 RF 
model on the cubic lattice. 

Here we present calculations for the K = 3 RA model on the Bethe lattice, partially 
motivated by the surprising result in the original simulations that no transition was found 
for this unrestricted model between the polymerized phases, while at least in the K = 2 
case the Bethe lattice calculations revealed no qualitative differences in the phase diagrams 
of the RA and RF models, both similar to the one found for the SASAW's model. In section 
HI] we define the model in more detail and present its solution on the Bethe lattice in terms 
of recursion relations. The thermodynamic behavior of the model is determined by the fixed 
points of the recursion relations, together with a bulk free energy which is useful to locate 
the coexistence loci, and these results may be found in section IIHI Final discussions and 
the conclusion are presented in section [TV] 



II. DEFINITION OF THE MODEL AND SOLUTION IN TERMS OF RECUR- 
SION RELATIONS 

n 

We study the MMS-RA model proposed by Krawczyk et al in [11] in the core of a Cayley 
tree with arbitrary coordination number q. In this model, self- and mutually avoiding walks 
are considered but the excluded volume condition is relaxed, so that each site of the tree may 
be occupied by up to K = 3 monomers, or, equivalently, each lattice site way be visited up 
to three times by the walks. No other restriction is imposed in the model, so that immediate 
reversals of the walk are allowed (RA model), differently from the more restrictive model 
studied in 6[ where immediate reversals are forbidden (RF model), and thus a subset of the 
configurations of the walks considered here was included. 

As usual, the endpoints of the walks are placed on the surface of the tree. Like in the 
original model fll|. a walk is described by the sequence of sites that it visits, so that the 
monomers placed on the same site are considered to be indistinguishable. A statistical weight 
Ui is associated to a site occupied by i monomers, with i — 1,2, 3. So, the grand-canonical 
partition function of the model will be given by: 

Y = ]T ufiupu,* (2) 

iVi,JV 3j JV 3 
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where the sum is over the configurations of the walks on the tree, while iV, , i = 1,2,3, is 
the number of sites visited i times by the walks. In Fig. [2J an example of a Cayley tree 
with three generations of sites is shown, as well as the contribution to the partition function 
which corresponds to the configuration of the walks in this case. 



FIG. 2: (Color online) A contribution to the partition function of the model on a Cayley tree with 



To solve the model on the Bethe lattice (the core of the Cayley tree) we start considering 
rooted subtrees, defining partial partition functions (ppf 's) for them, where we sum over all 
possible configurations of the chains for a fixed configuration of the root of the subtree (this 
is the reason for calling the partition functions partial). We thus define fourteen partial 
partition functions gi, i = 0, 1, ... , 14, shown in Fig. [3j The number of partial partition 
functions we need to define for the RA model is larger than the one we used for the RF model, 
where four root configurations were sufficient [6|, since more configurations are allowed in 
the present case. 

When we define the possible root configurations, with up to four polymer bonds on the 
root edge. It is important to notice that, since immediate reversals of the walk are now 
allowed, it is possible to have closed loops on the tree and the ppf's have to be carefully 
defined in order to avoid rings in the walks. This possibility does not exist in the RF model 





q = 4 and 3 generations. The weight of this contribution will be ti^o^wf. 
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on the Bethe lattice. When rings are allowed, even the universality class of a polymer 
model changes to a model with n — 1 components in the order parameter, that is, to the 
universality class of the Ising model. Therefore, in ppf 's with two or more bonds on the root 
edge, we need to distinguish between bond pairs that are connected (in earlier generations 
of the tree), or not. In Fig. [31 we have four ppf's with two bonds in the root site (g>2 to 
g$), for example. In bond pairs not connected by horizontal lines in our notation, such as 
in the root configurations for g^ and g§ (taken two by two) and one pair of g\o, the two 
walks will never meet on the tree at a site in earlier generations, or, in other words, if the 
walks are followed all the way to the surface of the tree, they end at different surface sites. 
When we draw horizontal lines connecting two bonds at the root edge, it means that they 
are connected to the same site in earlier generations and there are three possibilities for this: 
1) The two bonds are connected to the same monomer of the root site (one line); 2) Both 
walks are connected to the same site in some earlier generation and visit the same sites of 
the tree (two lines); and 3) Same as case 2, but the walks visit different sites of the tree 
(three lines). In the last case, the two bonds are distinguishable, because the sequence of 
visited sites is different if we begin in one or other bond, although the visited sites are the 
same. These definitions are applied for every ppf with at least two bonds in the root site, 
and leads to the rather large number of ppf's we need to define. 

OQOOOOOO 
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FIG. 3: Illustration of the rooted subtrees which correspond to the partial partition functions. The 
meaning of the horizontal lines between bonds on the root edge is discussed in the text. 



We then proceed obtaining recursion relations for the ppf's, by considering the operation 
of attaching q — 1 subtrees with a certain number of generations to a new root site and edge, 
thus building a subtree with an additional generation. Below the recursion relations are 
presented. In general, the partial partition function g[ with an additional generation is the 



S 



sum of contributions involving the parameters of the model and partial partition functions 
The primes denote the partial partition function of the subtree with one more generation. 
Whenever appropriate, the contributions to the sums begin with a product of two numerical 
factors, the first of which is the multiplicity of the configuration of the incoming bonds and 
the second is the multiplicity of the connections with the monomers located at the new root 
site. In the expressions below, /j = (^)go, where a = q — 1 is the ramification of the tree. 
The recursion relations for the 14 ppf 's are: 

9o = fo + ui[f 2 gl + fig 2 ] + W2L/4 x 3g* + 3/3 x 3glg 2 + Sf 3 gl(g 3 + g A + 2g 5 ) + 2/ 2 x 3 gi g 6 + 
2/2^1(^7 + 9s + 2^9) + h x 3#2 + 2/2^2(^3 + 94 + 2g 5 ) + /i#i ] + w 3 [/ 6 x 15gf + 
5/5 x Vzg{g 2 + 6/4 x 15gfg 2 + 12/ 4 x 6gfg 2 (g 3 + 94 + 2# 5 ) + 6/3 x l%i#2# 6 + 
6/3 x 6gMg 7 + g 8 + 2# 9 ) + 5/ 5 x 6g$(g 3 + g A + 2g 5 ) + 6/4 x 2g 2 1 (g 3 + g 4 + 2g 5 ) 2 + 
6/3 x 6 gi (g 3 + 94 + 2g 5 )g 6 + 6f 3 x 2 gi (g 3 + g 4 + 2^ 5 )(^ 7 + # 8 + 2# 9 ) + 4/ 4 x 15^# 6 + 
4/4 x 6^(37 + # 8 + 2^9) + 3/3 x 6^^io + 3/ 3 <?i 2 (<?ii + 2^i 2 + 4# 13 ) + 2f 29l g u + h x 15# 2 3 + 
3/3 x 6^(^ 3 + 94 + 2g 5 ) + 3/3 x 2g 2 (g 3 + 94 + 2g 5 ) 2 + 2/ 2 x 6# 2 £io + 2/2^2(^11 + 2#i 2 + 4# 13 ) + 
2/2 x 2((?3 + 54 + 2g 5 )g w + h x 15(? 6 2 + 2/2 x 6# 6 ((? 7 + g 8 + 2# 9 ) + / 2 x 2{g 7 + g 8 + 2g 9 ) 2 }; (3) 

g'i = uifigi + w 2 [/s x 3^ + 2/2 x 3#i£ 2 + 2f 2 gi(g 3 + #4 + 2# 5 ) + f\ x 3# 6 + /i(# 7 + #s + 2# 9 )] + 
w 3 [/ 5 x 15^ + 4/4 x 15^2 + 3/3 x 15^i^| + 6/3 x 6 gi g 2 (g 3 + g 4 + 2g 5 ) + 
4/4 x 6gl(g 3 + g 4 + 2g 5 ) + 3/ 3 x 2 9l (g 3 + g A + 2g 5 ) 2 + 3/ 3 x 15g 2 g 6 + 
3/3 x 6gl(g 7 + g 8 + 2g 9 ) + 2f 2 x 6 gi g l0 + 2f 2gi (g u + 2g 12 + Ag 13 ) + 2/ 2 x lhg 2 g 6 + 
2/2 x Qg 2 (g 7 + 9s + 2^9) + 2/2 x 6(g 3 + g A + 2g 5 )g 6 + 2f 2 x 2(g 3 + g 4 + 2g 5 )(g 7 + g s + 2g 9 ) + 

figu}; (4) 

9 2 = u 2 [f 2 gl + hg 2 ] + u 3 [U x 6gt + 3/ 3 x 6gfg 2 + 3/ 3 x 2g 2 1 (g 3 + g 4 + 2g 5 ) + 2f 2 x 6 gi g 6 + 
2/2 x 2 9l (g 7 + g 8 + 2g 9 ) + f 2 x Qg 2 2 + 2/ 2 x 2g 2 (g 3 + 94 + 2g 5 ) + f\ x 2(? 10 ]; (5) 
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wi + (Mh9i + hgt] + ^[fa x + 3/ 3 x 30102 + 3/301 (#3 + #4 + 20 5 ) + 2/ 2 x 30i0 6 + 
2/201 (07 + 0s + 20 9 ) + h x 30 2 2 + 2/202 (03 + 04 + 20 5 ) + /i0io] ; (6) 

= W2/l(03 + 04) + W 3 [3/ 3 0i(0 3 + 4 ) + 2/ 2 0i(0 7 + 8 ) + 2/202(03 + 04) + /i0io + /l01l(]?) 
5 = W 2 /l0 5 + W 3 [3/30105 + 2/ 2 0109 + 2/ 2 0205 + fl{gz + 04 + 20 5 ) 2 + /i(0i 2 + 20i 3 )]; (8) 



% = ^ 3 [/ 3 0i 3 + 2/20102 + /i0 6 ]; (9) 

W 2 /l01 + W 3 [/3 X 30 X 3 + 2/ 2 X 30102 + 2/20l(03 + 04 + 20 5 ) + /1 X 30 6 + /i(0 7 + 8 + 2#Jj) 

08 = W 3 [2/201 (03 + 04) + /l(07 + 08)]; (11) 

09 = u 3 [2/20105 + /i0g]; (12) 

0io = ^3[/ 2 0i + A02]; (13) 

g'u = u 2 + ^ 3 [/ 2 0i + A02]; (14) 

^2 = W 3 /l(0 3 +0 4 ); (15) 
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g'13 = ^3/1^5; 



(16) 



g'u = u 3 figi- (17) 

The partial partition functions are expected to grow exponentially with the number of 
iterations, so we define ratios of them, which usually remain finite in the thermodynamic 
limit. Furthermore, we notice in the above equations that some ppf's only appear in sums, 
they are (g 3 + # 4 + 2g 5 ), (g 7 + g 8 + 2g 9 ) and (g u + 2g 12 + 4# 13 ). Thus, it is convenient to 
define the following ratios: 

O 91 d 92 D (^3 + ^4 + 2^5) D #6 

fl'o go go go 

Rs= (97 + 91 + 29,) R9 = g 2 o ^ = ( 9 n + 2g 12 + Ag 13 ) ^ ^ = ^ 
fi'o fi'o fi'o fi'o 

From the recursion relations for the ppf's, similar expressions may be obtained for the 

ratios. Denoting the binomial coefficients as bi = (^) , the recursion relations for the ratios 

are: 



R[ = -[ou 1 b 1 R 1 +uj 2 (3b 3 Rl + 6b 2 R 1 R 2 + 2b 2 R 1 R 3 + 3b 1 R 4 + b 1 R 5 ) + 
u 3 (15b 5 Rl + mb 4 R\R 2 + A5b 3 RiR 2 2 + 3663^1^2^3 + 2Ab 4 R\R 3 + 
6b 3 RiRl + A5b 3 RjR A + 18b 3 RjR 5 + Ub^Re + 2b 2 R 1 R 7 + 30b 2 R 2 R A + 
12b 2 R 2 R 5 + 12b 2 R 3 R 4 + Ab 2 R 3 R b + b.Rs)}; (19) 



K = p[u 2 (b2Rl + b 1 R 2 )+uj 3 (6b 4 Rl + 18b 3 RlR 2 + 6b 3 R 2 1 R 3 + 12b 2 R 1 R 4 + 
4b 2 R x R 5 + 6b 2 R 2 2 + 4b 2 R 2 R 3 + 2b 1 R 6 )}; 



(20) 



K = + ^(hRj + b t R 2 + b t R 3 ) + uj 3 (3b 4 Rj + % 3 R\R 2 + Qb 3 R\R 3 + 

GhRiR* + IhRiRs + Zb 2 R\ + Ab 2 R 2 R 3 + 2b 2 R\ + 2b 1 R 6 + hRr)]] (21) 
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K = ^ M? + 2fe 2J R lJ R 2 + hR4; (22) 

i?s = -[u 2 b 1 R 1 +uj 3 (3b 3 R 3 l + 6b 2 R 1 R 2 + Ab 2 R 1 R 3 + 3b 1 R 4 + 2b 1 R 5 )}; (23) 

^6 = ^M? + Wk]; (24) 

R' 7 = -^[ou 2 +uj 3 (b 2 R 2 + b 1 R 2 + 2b 1 R 3 )}] (25) 

R' 8 = ^hR,. (26) 

The denominator D is defined as: 



D = 1 + wi(6 2 i?? + 6ii? 2 ) + w 2 (36 4j Ri + 963^1^2 + 363^1^3 + 662^4 + 2b 2 R 1 R 5 + 
3b 2 R 2 2 + 2b 2 R 2 R 3 + hR 6 ) + w 3 (15fc 6 #i + lhb 5 R\R 2 + 9064^2 + 72b 4 R\R 2 R 3 + 
mb 3 R x R 3 R 4 + 9063^1^2^4 + ZQhR^Rs + m b R\R 3 + ^i? 2 /^ + nb^R^ + 
mb 4 R\R 4 + 24b A RlR 5 + m 3 R\R & + 3b 3 RlR 7 + 2b 2 R 1 R 8 + 156 3 ^ 2 + 18b 3 R 2 2 R 3 + 
6b 3 R 2 R 2 3 + 12b 2 R 2 R 6 + 26 2 i? 2 i? 7 + 46 2j R 3 ^6 + 15& 2 ^ + 12b 2 R 4 R 5 + 2b 2 R\). (27) 

The grand-canonical partition function of the model on the Cayley tree may be obtained 
if we consider the operation of attaching q subtrees to the central site of the lattice, similar 
to the one used for deriving the recursion relations for the ppf's. The result will be: 

Y = g*(l + P + Q + S), (28) 
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where: 

P = ^(ciRj + Cl R 2 ); (29a) 
Q = u 2 {3c 4 B* + 9c 3 RfR 2 + 3c 3 i^# 3 + 602^1^ + 2c 2 R x R s + 

302^ + 202^3 + 01^6); (29b) 
S = lo 3 (15c 6 RI + 75c 5 RiR 2 + Wc^Rl + 12c±R\R 2 R 3 + 3Qc 3 RiR 3 R A + 
9003^1^2^4 + 3603^1^2^5 + 30c 5j R^ 3 + \2c A R\Rl + l2c 3 R 1 R 3 R 5 + 
mc A R\Ri + 2A Ci R\R b + 1803^1^6 + 3c 3 R\R 7 + 2c 2 R!R 8 + 15c 3 ^ + 
18c 3J R^3 + Qc 3 R 2 Rl + l2c 2 R 2 R 6 + 2c 2 R 2 R 7 + Ac 2 R 3 R 6 + \hc 2 R\ + 
12c 2 R A R 5 + 2c 2 R 2 5 ). (29c) 

where Cj = (^). We notice that the contributions to P, Q, S correspond to placing one, 
two and three monomers on the central site, respectively. Using the expressions above, we 
may obtain the densities at the central site of the tree, considering the configuration of this 
site for each contribution to the grand-canonical partition function |28j The density of sites 
occupied by one (pi), two (p 2 ) and three (p 3 ) monomers are given, respectively, by: 

Pl = l + P + Q + S ' (30a) 

P2 = i + p + g + s ' and (30b) 
p3 = i + p S +Q + s- (30c) 

The Bethe lattice solution of the model is defined by its thermodynamic behavior in the 
core of the tree, represented by the densities just defined. The total density of monomers 
on the Bethe lattice, that is, the total number of monomers divided by the number of sites, 
is p = pi + 2p 2 + 3p 3 , and will be in the range < p < 3. 



III. THERMODYNAMIC PROPERTIES OF THE MODEL 

The thermodynamic phases of the system on the Bethe lattice will be given by the stable 
fixed points of the recursion relations, which are reached after infinite iterations of the 
recursion relations and thus correspond to the thermodynamic limit. We find three different 
stable solutions for the fixed point equations R[ = Ri, associated to one non-polymerized 
phase (NP) and two polymerized ones (PI and P2). 
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The NP phase is characterized by the fixed point Ri = for all i, excluding % = 3 and 
i — 7. These last two may be obtained solving the equations: 



2b 2 co 3 R 2 3 + (huz + 2b\ul - 1)R 3 + (u! + &iw 2 w 3 ) = 0, (31) 

and 

R 7 = u 2 + 26iw 3 i? 3 . (32) 

The quadratic equation for R% can be easily solved, but the explicit expression for R^ p 
is too large to be shown here. Looking the equations [29] and [30] we see that p = in the NP 
phase, as expected. 

In the two polymerized phases all ratios are non- vanishing and, in order to obtain the fixed 
point values, we have to iterate the recursion relations or solve the fixed point equations 
numerically. It is important to remark that the metastable phases with double e triple 
occupation of sites, that appear in the RF model [6(, are absent here. As is discussed in 
[l3| , the immediate reversal of the walk makes the probabilities of find these configurations 
in the central site vanish. 

The stability limits of all phases are obtained calculating the jacobian of the recursion 
relations: 

'«=(§)■ < 33 > 

A fixed point is stable if the dominant eigenvalue of the jacobian has a modulus smaller 
than unity, and the stability limit of the corresponding phase (spinodal) is located in the 
loci where this modulus becomes equal to unity. 

In order to find the coexistence surfaces in the phase diagrams, we obtain the free en- 
ergies of the thermodynamic phases of the model. This free energy may not be calculated 
directly from the partition function Eq. (|28|) . since it refers to the whole Cayley tree and, 
remembering that in the thermodynamic limit the number of surface sites correspond to a 
finite fraction of the total number of sites, reflects the influence of the surface. The free 
energy per site in the core of the tree, which corresponds to the Bethe lattice solution, may 



be calculated following the Gujrati's argument 16j, which was also derived in another way 



in jo]. The result for the grand-canonical free energy per site on the Bethe lattice (divided 
by ksT) is: 

(j> b = -~ [q In (D) - (q - 2) In (1 + P + Q + S)} (34) 
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Using the spinodals to find the continuous transitions and the free energy to determine the 
coexistence surfaces we built the whole phase diagram of the system. Before presenting the 
complete three-dimensional phase diagram, in the space defined by the statistical weights 
Ui, 0U2, and U3, it is instructive look at the thermodynamic behavior in some cuts of the 
parameter space. All results presented below are for a lattice of coordination q = 4 and 
qualitatively identical 1 
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FIG. 4: (Color online) Phase diagram for CJ3 = 0. The red curve (above the tricritical point-TCP) 
is a first order transition and the black line (below the TCP) is a continuous transition between 
the NP and P2 phases. The TCP is indicated by the green square. 



The diagram for U3 = (K = 2) is shown in Fig. HI For small values of 0J2, we 
find a continuous transition, between the NP e P2 phases, which ends up at a tricritical 
point (TCP) located at ^^^=2) = Q mj and jtcp{k=z) = QAUG ( for q = A y Above 

the tricritical point the transition becomes discontinuous. It is important to stress that this 
particular case (co 3 = 0) was studied by one of the authors in [13| , considering distinguishable 
monomers and using the natural initial conditions method to find the coexistence lines. 
There, only a discontinuous transition was found between the non-polymerized phase and a 
polymerized one (called P2 here). We notice that the diagram of 13] changes if Gujrati's 
prescription is used to obtain the bulk free energy and the coexistence lines are determined 
using this free energy. Since this latter procedure has a better fundamentation than the 
earlier based on natural initial conditions [6], the results provided by it are more reliable. 
We also notice that the phase diagram found here is very similar to one obtained for the RF 



model in |6j. However, the RF tricritical point was located at uif =1/3 and uj. 



,RF 



1/9 
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6), which is far from the location found here. A generalization of the model for K = 2, with 
the RA and RF models as particular cases, shows a line of tricritical points joining the two 
ones for both models Ijj]. 

For increasing values of U3, the qualitative behavior of the phase diagrams, in the (ui, 002) 
plane, is similar to the one depicted in Fig. SJ Thus, there is a tricritical point (TCP) line 
in the three dimensional phase diagram. This line ends up at a multicritical point (MCP), 
located at (cof WP = 0.2700819, CP = 0.0182769, CP = 0.0450771), so the TCP line 

MCP -n^ , , ^ , MCP 



J 3 

lies in the region < U3 < u^ 1VP . For U3 > uj^ 1up ', more complex diagrams are found, 



as will be discussed below. The location of this multicritical point may be determined by 
noting that it corresponds to a higher order NP root of the fixed point equations. This will 
be discussed in some d 
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FIG. 5: (Color online) Phase diagram for 0J2 = 0. At the red (left of the critical endpoint-CEP) and 
the magenta (right of the CEP) curves we find first order transitions between the polymerized phase 
(PI) and the non-polymerized phase (NP) and between both polymerized phases, respectively. 
The black line (below the CEP) corresponds to continuous transitions between the NP and the 
polymerized phases. The blue (square) and magenta (circle) points are the critical endpoint (CEP) 
and a critical point (CP), respectively. All these features are better seen in the inset, which shows a 
small region close to the point where the discontinuous (NP-P2) and continuous (NP-polymerized) 
transitions lines meet. 



A rich phase diagram was found in the U2 = plane, as show in Fig. [5j For small values 
of uj\ there is a first order transition between the NP em P2 phases, ending at a critical 
end-point (CEP), located at oof EP « 0.2823 and oj$ ep « 0.0518. In a tiny region of the 
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phase diagram, where u>\ > Ui EP and CU3 



< u CEP 



we found the second polymerized phase 



(PI). The two polymerized phases (PI and P2) coexists in line which limits this region until 
a critical point is reached (at wf p ~ 0.2831 and oj^ p ~ 0.0516). This is shown in the inset of 
figure EJ Below the coexistence line (for u\ > uf EP and W3 < u^ EP ), there is a continuous 
transition line befrw 




0.050 



FIG. 6: (Color online) Densities as a function of W3, for ui\ = 0.2799 and L02 = 0.004. In the inset 
we show a detail of the region with discontinuous transitions. 

Comparing the phase diagram of Fig. [5] with the one of the RF model (for u 2 = 0) |(|, 
we can see that the qualitative picture is the same. However, in the RF model, the P1-P2 

n 

coexistence region is larger than the one for the RA model. Besides, as shown in [6J, for 
the RF model, the phase PI is characterized by p\ ^> p 2 ,P3 at the coexistence with the P2 
phase, where pi ~ p 2 ~ Pz- On the other hand, here we found that in the PI phase all 
densities are of the same order as in P2, and thus the two phases have the same symmetries. 
In figure El we show the densities (defined in Eq. [SU]) for increasing values of u> 3 , with U\ 
and u 2 fixed in the P1-P2 coexistence region. In face of this result, we can conclude that 
the restriction imposed in the RF model changes the nature of the PI phase, which became 
approximately a SAW in that case, with dominance of sites with a single monomer, while 
when immediate reversals (RA model) for the walks are allowed, the PI phase behaves like 
a regular polymerized phase in the MMS model. We advance that this difference of the PI 
phase in the two models may explain the different regimes found in the phase boundaries in 



the canonical simulations of Krawczyk et al 
below. 



111 ]. This point will be discussed in more detail 
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The (001,003) phase diagrams are similar to the one in Fig. [5j for all 002 smaller then the 
multicritical point value (002 < oo^ 10 )- Therefore, we have lines of CP and CEP in the whole 
phase diagram and a P1-P2 coexistence surface between these lines. Thus, there exist a 
NP-P2 coexistence surface and a NP-polymerized critical surface in the three dimensional 
diagram. For increasing values of 002, the CP line gets closer to the CEP line, making the 
numerical determination of their locations very hard. At the multicritical point these lines 
meet, together with the TCP line. When oo 2 MCP < 002 < oo^ CP 2 ^ ', the tri critical point line 
crosses the (001,003) plane and the oo\ x 003 diagrams resembles the one shown in Fig. HJ 

In Fig. we show several diagrams, in the (002,003) plane, for different fixed values 
of oo\. For 001 = (Fig. [TJa)), there is only a NP-P2 coexistence line and, for oo\ < 
^tcp(k-2) ^ g j m jj ar cii a g ram s are obtained, forming the NP-P2 coexistence surface. In the 
range < Wi < oof ICP , the NP-P2 transition may be continuous or discontinuous, 

and both transition lines meet at a tri critical point. In Fig. [7J (b) we show an example of 
these diagrams for oo\ = 0.23. Finally, for oo\ > oo^ CP , we find the same behavior of the 
(oo\, 003) diagram for small 002 (see Fig. [5]), with two coexistence lines (NP-P2 and P1-P2) 
which meet at a critical end-point, where the line of continuous transition between the NP 
and the polymerized phases ends (see Fig. (c)). 

Again, the diagrams found here with fixed oo\ are similar to those of the RF model jg]. 
The main difference is that the tricritical and critical end-point lines of the RA model are 
functions of the three parameters (oo\, 002 and 003), while in the RF model these lines lie in 
the plane oo\ = 1/3. In the same way, in the RA model the NP-polymerized continuous 
transition appear as a curved surface, while in the RF model, it is located in the plane 
u x = 1/3. 

A sketch of the whole three-dimensional phase diagram is shown in Fig. [HI summarizing 
the features we have discussed above. Like discussed above, the CEP and CP lines are very 
close in the phase diagram and we can not see the two lines in the figure Therefore, we 
show only the CEP line in the diagram, but it is important to keep in mind that there is 
also a CP line in the neighborhood of this line. In particular, due to the existence of this 
additional coexistence surface between both polymerized phases, the NP-P critical surface 
and the NP-P coexistence surface do not meet tangentially at the CEP line, and the angle 
between the normal vectors to both surfaces at this line becomes larger as further we are 
from the multicritical point, where the tricritical, critical endpoint and P1-P2 critical line 
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FIG. 7: (Color online) Phase diagrams for: a) uj\ = 0.0, b) u)\ = 0.23 and c) u% = 0.275. The red 
and magenta (dashed) curves are NP-P2 and P1-P2 discontinuous transitions, respectively. The 
black line is the NP-polymerized continuous transition. 



meet tangentially. Therefore, this angle is largest when U2 = 0, as may be seen in Fig. |5j 
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FIG. 8: (Color online) Sketch of the three-dimensional phase diagram. The NP-P2 discontinuous 
transition surface (red, dashed lines) and the NP-polymerized critical surface (black, continuous 
line) are shown. The P1-P2 coexistence surface and critical line close to the CEP line are to small 
to be represented here. 

IV. FINAL DISCUSSIONS AND CONCLUSIONS 



Although the Bethe lattice solution of the RA model is close to the one of the restrictive 
RF model, the polymerized phases PI are very different in the two models. Moreover, in the 
RF model the continuous transition surface exists between the NP and PI phases only, while 
here the main part of this surface are between the NP and a polymerized phase that can 
not be identified as PI or P2, because they are above the critical point line. These results 
may explain the difference found in the canonical simulations of Krawczyk et al TjJ between 
the RA and RF models. For the RF model, Krawczyk et al [n| suggest a canonical phase 
diagram with discontinuous and continuous transition lines which meet at a tricritical point, 
between a SAW-like phase (sites occupied mainly by a single monomer) and a collapsed one 
(sites predominantly with two or three monomers). The location of the TCP is not denned 
precisely by the simulations, but the authors suggest that it is in the region of attractive 
interaction between monomers, namely, the first quadrant in the (Pi, P2) parameter space. 
The Bethe lattice solution of this model p shows that, in fact, there is a SAW-like phase (the 
critical surface NP-P1) and a collapsed phase (the coexistence surface NP-P2). However, the 
discontinuous and continuous transitions lines, suggested in the simulations, are a tricritical 
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and a CEP line in this approximation, respectively. On the other hand, no SAW-collapsed 



transition was found in the simulations of the RA model in [11] . This is in agreement 
with our results, because here the critical surface NP-polymerized does not lead to a SAW- 
like phase in the canonical diagram. In contrast with the RF model, where the sites are 
predominantly visited by one monomer, here the densities in phase PI depends only on the 
statistical weights like in phase P2. Thus, both the critical surface NP-polymerized 
and the coexistence surface NP-P2, are associated with collapsed phases, since the sites are 
occupied predominantly by more than one monomer in both cases. The former leads to 
a collapsed phase with low density (CLD) of monomers and the last has a larger density 
(CHD). We believe that the similarity between these phases makes it difficult to distinguish 
them in the simulations, which could have lead to the conclusion of no transition for the 
canonical RA model in [11] . 



In order to compare our grand-canonical resu 



ts for the RA with the canonical ones 



obtained in the simulations for the RF model 11], we map our grand-canonical diagram 
into canonical one. As was discussed in the canonical variables used in the simulations 
11] are related to the Boltzmann weights of our solution as: 



01 = In 



u 2 



and 



02 = In 



^3 



(35) 



and in the canonical formalism we are always restricted to the boundary of the NP phase. 
The canonical phase diagram which we found in the 0\ , 02 parameter space is shown in Fig. 
[9j Like discussed above, we found two collapsed phases with high (CHD) and low (CLD) 
monomer densities that are related to the grand-canonical NP-P2 coexistence surface and 
NP-polymerized critical surface, respectively. 

The CLD-CHD transitions are always continuous, but of different types: for values of 02 
below the multicritical point there is a tricritical line and above this point a critical end- 
point line separates the two phases. The multicritical point is located at 0\ = —1.3840554 
and Bo = 0.8277081. The same behavior was found in the Bethe lattice solution of the RF 

n 

model [6], but there the whole TCP line lies in the region of negative values for 2 negative 
axis (0i = 0, 02 < 0) and the MCP is placed at the origin (0i —02 — 0). Curiously, here 
the MCP is located in a region with repulsive interaction between monomers at same site, 
when only two monomers are present. Thus, in this region, sites occupied by two monomers 
are penalized, while sites with one or three monomers are favored. In fact, in the Fig. |6l 



21 



1.0- 



MCP 



P 2 0.4 



-0.2- 



0.0- 



0.2- 



0.8- 



0.6- 



NP-Polymerized 
(CLD) 




NP-P2 " 
(CHD) '_ 



-3 



-2.5 



-2 



-1.5 



-1 



-0.5 







0.5 



FIG. 9: (Color online) Canonical phase diagram. The green curve (below the MCP, represented 
by the black circle) is the tricritical line and the blue (above the MCP) is the CEP line. On the 
line p 2 = 3/3i (dashed) the MMS-RA model is related to the SASAW's model. 

that show the densities in a region close to the MCP, we can see that pi, ^> p2- 

A connection can be established between the MMS and the SASAW's models, as was 
already discussed in If we suppose that only monomers located at the same site have 
through an attractive pairwise interaction of energy — e, in the grand-canonical ensemble 
we should have u\ = z, w 2 = z 2 u, and o>3 = 3z 3 u, where, as before, uj = exp(/3e). In the 
canonical situation, from Eqs. ( 1351) . this leads to 02 — 3/?i = 3a;. This line is shown in Fig. 
IHJ and it crosses the tricritical line, so that the collapse transition in the subspace of the 
parameter space where the MMS-RA model on the Bethe lattice with K = 3 is related to 
the SASAW's model is tricritical, as it is also in the SASAW's model. 

In conclusion, to study the thermodynamic behavior of models for complex fluids such 
as the one considered here, for which exact solutions are usually difficult to obtain, we 
think it is useful to combine numerical simulations with approximate analytic solutions. In 
particular, the findings in this work suggest that the qualitative behavior of the MMS model 
without restrictions (RA) may be similar to tho one of the restricted model (RF). Also, on 
the Bethe lattice, the MMS model shows a behavior which is close to the one of the standard 
SAW's model for the collapse transition of polymers. Of course one has to be aware that 
Bethe lattice solutions, as all mean-field like approximations, overestimate the interactions, 
and therefore may predict phase transitions in situations where better approximations or 
exact solutions find none, but in our opinion the results found on the Bethe lattice for the 
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MMS model suggest that additional investigations using simulations or other more precise 
techniques are desirable. 
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Appendix A: Determination of the location of the multicritical point 

The multicritical point may be defined as the common point of the lines of tricritical 
points and critical endpoints, as may be seen in the full phase diagram depicted in Fig. 
El This definition, however, does not lead directly to a precise algorithm to determine the 
location of the multicritical point, due to the rather intricate nonlinear fixed point equations 
which need to be handled for this purpose. Sometimes, in Bethe or Husimi lattice solutions, 
it may be possible to reduce the fixed point equations to a polynomial, and then the higher 
order transition loci can be identified with the higher order roots of the NP fixed point, and 
example of this procedure is described for the particular case K = 2 of the RF-RA model 
in 14j. We were not able to pursue accomplish this calculation in thepresent case. 

Another possibility, which was used for the K = 3 RF model in [6[, is to assume that, 
close to the NP fixed point, the remaining ratios may be expanded in powers of one of them. 
By expanding the fixed point equations in powers of the chosen ratio, one then requires the 
higher order transition loci to be a higher order root of the resulting set of equations in the 
parameters of the model and the expansion coefficients. In the present case, we expanded 
the remaining ratios in powers of Ri and found a consistent solution of this kind by requiring 
the multicritical point to satisfy the expanded fixed point equations up to order Rf. This 
rather high order is necessary due to parity effects we found in the expansions of the ratios. 
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We assumed that, up to order Rf, the ratios may be expanded as follows: 



R 2 = a 22 R\ + a24-Ri, (Ala) 

#3 = a 3 o + «32^i + a 34 Rj, (Alb) 

Ra = a 43 R\ + a Ab R\, (Ale) 

R 5 = a 51 R! + a 53 Rl + a 55 #i, (Aid) 

R 6 = a G2 R\ + «64-Ri, (Ale) 

R7 = a-70 + a-j 2 R\ + o-iaRi, (Alf) 

R 8 = a 81 Ri + a 83 Rl + a 8S R\. (Alg) 



Now these expressions are substituted into the 8 fixed point equations which are obtained 
by imposing R[ = Ri in the recursion relations for the ratios Eqs. f lT9H26|) . Expanding these 
fixed point equations up to order R\ , using a algebra software for this task, we obtain several 
expansion coefficients of the fixed point equations. Considering the parity of the expansions 
shown in Eqs. ( 1A1I) and denoting by the coefficient of order R\ in the fixed point 
equation obtained from the recursion relation for we are lead to 21 equations Cjj = 0, 
corresponding the the coefficients (1,1), (1,3), (1,5), (2,2), (2,4), (3,0), (3,2), (3,4), (4,3), 
(4,5), (5,1), (5,3), (5,5), (6,1), (6,3), (6,5), (7,0), (7,2), (7,4), (8,1), (8,3) and (8,5). In 
particular, C 3i0 = and CV = lead to the pair of equations (13 1|I32 j) for the fixed point 
values of the NP phase presented above. The complete set of equations is too long to be given 
here. Solving this set of nonlinear algebraic equations for the 18 expansion coefficients and 
the 3 statistical weights Ui, we obtain the result: a 22 = 0.4882335108, a 2 4 = 0.2869523407, 
a 30 = 0.3329767741, a 32 = -0.03809592295, a 34 = 0.1519415331, a 43 = 0.2048247034, 
a 45 = -0.2633449211, a 51 = 0.3220483966, a 53 = 0.1813751089, a 55 = -0.6519478195, 
a 62 = 0.2012557404, a 64 = -0.2612169532, a 70 = 0.1083347046, a 72 = 0.03872548434, a 74 = 
-0.1169008783, a 8 i = 0.1352312919, a 83 = -0.2015959772, and a 85 = 0.1269313668 for the 
expansion coefficients and ui = 0.2700819945, u 2 = 0.01827694593, and u 3 = 0.04507709731 
for the statistical weights. We did some numerical calculations for the fixed point values of 
the ratios in the neighborhood of the multicritical point, and found that they are consistent 
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with the expansion coefficients we found. 
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